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Abstract 

We analyze a recent application of homotopy perturbation method to some heat -like 
and wave-like models and show that its main results are merely the Taylor expan- 
sions of exponential and hyperbolic functions. Besides, the authors require more 
boundary conditions than those already necessary for the solution of the problem 
by means of power series. 



1 Introduction 



In the last years there has been great interest in the power-series solution 
of physical problems by means of homotopy perturbation method [1-5]. In a 
series of papers I have shown that the homotopy perturbation method, the 
Adomian decomposition method and the variational iteration method have 
produced the poorest and most laughable papers ever published [6-12]. Notice, 
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for example the application of homotopy perturbation method to obtain the 
Taylor expansion of exponential functions of the form e*"* [4] . Curiously, some 
journals believe that such contributions may be of interest for the scientific 
community. 

The purpose of this article is to discuss a recent paper by Ozi§ and Agirseven [5] 
who have recently applied the homotopy perturbation method to heat-like 
and wave-like equations. We compare their approach with the straightforward 
power-series one. 



2 Heat like models 



Ozi§ and Agirseven studied particular cases of the equation [5] 
du{r, t) 



dt 



Luir,t) + f{r) 



where L is a linear operator and /(r) a differentiable function. As a first 
approximation we can naively try a time-power-series solution of the form 

oo 

«(r,t)=^w,(r)t^- (2) 

If we substitute equation ([2]) into equation ([1]) we obtain a recurrence relation 
for the coefficients of the expansion ([2]): 

Un+i = -^—r (LUn + f6no) ,n = 0,l,... (3) 

where bij is the Kronecker delta. Our first observation is that we only need 
the coefficient 

Uo(r) = ^/(r,0) (4) 



in order to obtain the remaining ones by means of the recurrence relation ([3]). 
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The first example studied by Ozi§ and Agirseven [5] is a particular case of 
equation ([T]) with 

If we apply the recurrence relation ([3]) with ^0(2;) = u{x,0) = we obtain 
the coefficients of the time-power series for the exact solution 

u{x,t)=x^e^ (6) 



that one easily obtains by the method of separation of variables. It is amaz- 
ing that Ozi§ and Agirseven [5] require the additional boundary conditions 
M(0,t) = and u{\^^) = e*. On the other hand, the power-series method 
reveals that the initial condition completely determines the solution obtained 
by those authors and that the remaining boundary conditions are redundant. 
But the most important fact is that those authors were allowed to publish an 
elaborated perturbation method for obtaining the Taylor series of e*!!! [5]. 

The reader may easily verify that the straightforward application of the re- 
currence relation enables one to solve all the other heat-like equations 
chosen by Ozi§ and Agirseven [5], and in all the cases one needs only the ini- 
tial condition (jl]) in order to obtain the solutions. On the other hand, Ozi§ and 
Agirseven [5] require some additional boundary conditions. For completeness 
we outline the results in what follows: 

Example 2: 

The authors require four Neumann boundary conditions and the initial one 
u{x,y,0) = y"^ [5]. The power series method completely determines the exact 
solution u{x, y, t) = cosh t + x"^ sinh t from only the initial condition. Notice 
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that in this case the authors were able to obtain the Taylor series for sinht 
and cosht which are far more difficult than the preceding feat. 

Example 3: 

The authors require six Neumann boundary conditions and the initial one 
u{x,y,z,0) = [5]. The time-power series completely determined by the 
initial condition clearly converges towards the exact solution u{x,y,z,t) = 
x^y'^z^{e^ — 1). In this case the authors clearly show that the homotopy per- 
turbation method is able to provide the Taylor series for e* — 1 about t = 0. 



3 Wave-like models 

Ozi§ and Agirseven [5] also studied some wave-like equations that are partic- 
ular cases of 

^^ = Ln(r,t) + /(r) (9) 

If we substitute the time-power series ([2]) we realize that one can obtain the 
coefficients from the recurrence relation 



Un+2 = -1— 

(n 



(Lun + fSno) , n = 0, 1, . . . (10) 
1 n + 2 m ^ 



In this case if we have the first two coefficients 

du{r, t) 



Uo{r) = u{r,0), Ui{r) 



dt 



(11) 



we can obtain the remaining ones. The example 4 studied by Ozi§ and Agirseven 
[5] is a particular case of equation with the linear operator and function 
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given by equation ([5]). If we choose the coefficients uo{x) = x and ui{x) = x"^ 
proposed by Ozi§ and Agirseven [5] we easily obtain the remaining coefficients 
of the time-power series of the exact solution 

u{x,t) = X + x^ svah.{t) (12) 

Once again we see that we do not need the other boundary conditions m(0, t) = 
0, and u{l,t) = 1 + sinh(t) already required by Ozi§ and Agirseven [5]. In this 
case the authors were able to produce the Taylor series for sinh(t). 

We draw the same conclusion regarding the other wave-like models studied 
by Ozi§ and Agirseven [5] as shown in what follows: 

Example 5: 

The authors require four Neumann conditions in addition to the initial ones 
Uq{x, y) = x^, and Ui{x, y) = y'^ [5]. One easily verifies that the straightforward 
power-series method yields the exact solution u{x, y, t) = x"^ cosht + y^sinht 
from only the two initial conditions already indicated. Once again, the authors 
were able to obtain the Taylor series for cosht and sinht. 

Example 6: 

The authors made use of six boundary conditions and the two initial ones 
Uo{x, y, z) = 0, and ui{x, y,z) = x"^ + y'^ — [5]. In this case we verify that the 
power-series method with the two initial conditions completely determine the 
exact solution it(x, ?y, z, t) = (x^ + y'^){e^ ~ 1) + -2^(e~* — 1)- As the reader may 
have already guessed, Ozi§ and Agirseven [5] were able to derive the Taylor 
expansions for e* — 1 and e"* — 1 about t = 0. 
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4 Conclusions 



We have seen that the homotopy perturbation method proposed by Ozi§ and 
Agirseven [5] produce the Taylor series about t = of the solutions of some 
differential equations. Curiously, those authors require more boundary condi- 
tions than the ones already necessary for the straightforward application of the 
power-series approach. The main results in their paper are Taylor expansions 
of exponential and hyperbolic functions of time about t = 0. It is amazing that 
a supposedly respectable journal accepts such contribution. However, I may 
be mistaken and in a near future all the courses on calculus will be teaching 
homotopy perturbation theory instead of Taylor series. 
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